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We make use of recent results in effective theory and higher-order perturbative calculations to 
improve the theoretical predictions of the top-quark pair production cross section at hadron colliders. 
In particular, we supplement the fixed-order NLO calculation with higher-order corrections from soft 
gluon resummation at NNLL accuracy. Uncertainties due to power corrections to the soft limit are 
estimated by combining results from single-particle inclusive and pair invariant-mass kinematics. 
We present our predictions as functions of the top-quark mass in both the pole scheme and the MS 
scheme. We also discuss the merits of using threshold masses as an alternative, and calculate the 
cross section with the top-quark mass defined in the IS scheme as an illustrative example. 



PACS numbers: 14.65.Ha, 12.38.Cy 



I. INTRODUCTION 



The total tt production cross section is an important observable at hadron colliders such as the Tevatron and LHC. 
For instance, it provides information about the top-quark mass, which is an input for electroweak fits [2| used to 
constrain the mass of the Higgs boson. Extractions of the top-quark mass from the production cross section have 
the advantage that the perturbative calculations used in the analysis are carried out in a well-defined renormalization 
scheme for the top-quark mass. The pole mass as well as the MS mass have already been extracted from the production 
cross section at the Tevatron 0]. The use of a short-distance mass such as the MS mass is theoretically favored over 
thepole mass, which can be defined only up to a renormalon ambiguity of order Aqcd- Moreover, it has been proposed 
in [3| that the apparent convergence and scale uncertainties of the perturbative series for the total cross section is 
improved in the MS scheme compared to the pole scheme, already at low orders in the perturbative expansion. 

According to the QCD factorization theorem [J], the hadronic cross section can be obtained from the partonic one 
after convolution with parton distribution functions (PDFs). For tt production, it is conventional to express this 
factorization theorem in the form 



a{s,mt) = 2_. y ^ — -^y ( -'^/ ) ( ^'^/ ) ' (1) 



s 



where /i/ is the factorization scale, i,j G {q, q, g}, and the parton luminosities are defined as 

/■I dx 

ffij{y,l^f)= — h/N^{x,lif) fj/N^v/x^l^l) ■ (2) 

Jy X 

Here and throughout this letter the running coupling is defined in the MS scheme with five active (massless) flavors. 
The scaling functions fij are proportional to the partonic cross sections and can be expanded in powers of as/ir. 
Numerical results for the next-to- leading order (NLO) term have been known for over two decades [ll-Q, and more 
recently fully analytic results were obtained in Q. Predictions based on NLO calculations typically exhibit scale 
uncertainties larger than 10%. In order to further reduce the theoretical uncertainties to match the experimental 
precision, it is necessary to go beyond NLO. Consequently, the calculation of the next-to-next-to-leading order (NNLO) 
corrections is an active area of research. Many efforts have been made in the calculation of two-loop virtual corrections 
, one-loop interference terms [Tsl-flTj , and double real emission [T8l - [2l| . Due to very recent progress in developing 
a new subtraction scheme for double real emission in the presence of massive particles 22,-26.] , a calculation of the 
NNLO corrections now seems feasible, and its completion will be a major accomplishment. 

An important way to improve on the fixed-order calculations is to supplement them with threshold resummation 
[27, 28]. In such an approach, which in general works at the level of differential cross sections, one identifies a partonic 
threshold parameter which vanishes in the limit where extra real radiation is soft, and sums a certain tower of 
logarithmic corrections in this parameter to all orders in the strong coupling. Such resummed formulas neglect power 
corrections which vanish when the threshold parameter goes to zero, but these can be taken into account, up to a given 
order in the strong coupling, by matching with the fixed-order results. In this way, one obtains predictions which not 
only make full use of the fixed-order calculations, but also resum a class of logarithmic corrections to all orders. In 
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the limit where the higher-order corrections are dominated by these logarithmic terms, such a resummation is clearly 
an improvement. For this reason, many different implementations of threshold resummation have been considered in 
the literature, where the current frontier is next-to-next-to- leading logarithmic (NNLL) accuracy jsl. [29l - [35j . 

The purpose of this letter is to consolidate results for the total cross section based on threshold resummation in 
effective field theory. These are obtained from two different threshold limits for the differential cross section. The first, 
referred to as pair invariant-mass (PIM) kinematics, uses the top-pair invariant-mass distribution as the fundamental 
observable [s^, H^- The second, referred to as single-particle inclusive (IPI) kinematics, works at the level of the 



transverse-momentum or rapidity distributions of the top quark |35j . The total partonic cross section is obtained by 
integrating the resummed distributions over the appropriate phase space. In both kinematics, the predictions in the 
effective field-theory framework include resummation effects to NNLL order, and are matched with the fixed-order 
results at NLO in order to achieve NLO+NNLL accuracy for the total cross section. Such formulas can be evaluated 
numerically using specific values of the matching scales appearing in the effective-theory calculations, or otherwise 
re-expanded in a fixed-order expansion, in the form of approximate NNLO predictions. We have considered both 
scenarios in [s^ls^, finding that when the results in PIM and IPI kinematics are combined as in the present work, 
the numerical differences between the NLO+NNLL and approximate NNLO results are rather small. In this letter 
we focus on the approximate NNLO results for concreteness. 

In what follows, we briefly review the formalism, and then give our best numerical predictions for the total cross 
section in the pole scheme for the top-quark mass, in the form of numerical fits as functions of mt- We then explain 
how to convert the results to the MS scheme and present numerical results as a function of the MS mass. The results 
in both schemes are calculated with a Fortran implementation of our approximate NNLO formulas, which we include 
with the electronic submission of this letter. Contrary to Q, we do not find a strong improvement of the perturbative 
series in the MS scheme compared to the pole scheme. We suggest that the group of short-distance masses known as 
threshold masses (see e.g. |36|) may actually be the more appropriate choice, and we explicitly consider the case of 
the mass defined in the IS scheme (37j . However, we believe that the poor large-order behavior of the perturbative 
series in the pole scheme is unlikely to be of practical importance in the foreseeable future. 



II. INGREDIENTS OF THE CALCULATION 



Our calculation is based on [33|, |35|, where threshold resummation using soft-collinear effective theory (SCET) 
[ssl - liol l is applied to tt production at hadron colliders. The partonic scattering process is 



iiPi) + ]{P2) ^ i(P3) + t{pi) + X{k) 



(3) 



where j, j indicate the incoming partons and X is a partonic final state. In [33|, the invariant-mass distribution 
da/dM of the tt pair was considered, where = (pa + PaY , and the threshold region is defined by 2: = M'^ / s — > 1 
with s = {pi +^2)^ (so-called PIM kinematics). In [35|, on the other hand, the transverse-momentum and rapidity 



distributions of the top-quark were considered. In this case, the threshold region is defined by S4 = {p4 + k) 



- mt 







(IPI kinematics). In both the PIM and IPI threshold limits, the energy of the extra emitted gluons vanishes, and the 
partonic differential cross sections are dominated by singular distributions of the form 



In^e 



0, ...,2n- 1 



(4) 



where ^ = (1 — z)/y^ for PIM and ^ — S4,/{mt\/rn^~+T2) for IPI kinematics. 

In the threshold limit ^ — >■ 0, the partonic differential cross section can be factorized into a product of matrix-valued 
hard and soft functions, 



da{£,, fif) oc Tr 



(5) 



where we have suppressed the dependence on the other kinematic variables, and £'«($) is the energy of the soft gluon 
radiation, which is given by — in the partonic center-of-mass frame for PIM kinematics, and Es — mt ^ in the 
t rest frame for IPI kinematics. The hard function H is related to virtual corrections and is thus evaluated at ^ = 0. 
The ^ dependence of the cross section is encoded in the soft function S via the ratio Es{£,)/fJ,f. One can therefore use 
the renormalization-group (RG) equation of the soft function to resum the singular distributions in ^ to all orders in 
as ■ A technical complication is that the RG equation of the soft function is non-local. We solve this equation using 
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the technique of Laplace transforms [4l|, and then carry out the inverse Laplace transform analytically to obtain 
expressions in momentum space. The final result for the resummed cross section reads 



da{i,Hf ) oc exp[4a^*(/is, ^/)] Tr U{fih,fJ-s) H{fih) U'' {p.h, fJ-s) s{d,,,as{ps)) 



r(2?7) 



(6) 



where s is the Laplace transform of the soft function, U is an evolution matrix which resums large logarithms 
between the hard and soft scales fih and fis, and 77 = 2a^{^J■s, l^f) and a^4,{^s, l^f) are related to certain integrals over 
anomalous-dimension functions. Explicit expressions for these objects in RG-improved perturbation theory can be 
found in [s^. For NNLL accuracy, one needs the anomalous-dimension matrices for the hard and soft functions to 
two-loop order (computed in [42*1), as well as the hard and soft functions to one- loop order (computed in [ssl [35|). 
One must also specify a procedure for choosing the hard and soft scales. In the SCET approach, these are chosen 
such that the contributions of the NLO matching functions to the hadronic cross section are minimized; details can 
be found in [H, HI] . 

The effective theory is a powerful tool for separating physics at the formally very different scales /i^ and /is and for 
summing logarithms of their ratio. In practice, however, in areas of phase space where the differential cross section 
is largest, the numerical hierarchy between the two scales is only moderate. Therefore, an alternative to using the 
NLO-I-NNLL formulas directly is to re-expand them in a fixed-order series in as(/i/), constructing what are referred 
to as approximate NNLO formulas. The full NNLO correction to the cross section takes the form 



ln'"t 



Co 6(0 + m) 



(7) 



The approximate NNLO formulas derived from the NLO+NNLL results contain an unambiguous answer for the Dm 
coefficients in the limit ^ — 0. Only parts of the coefficient Co are determined, and the regular piece R is subleading in 
the threshold limit and not determined. There are thus ambiguities in what to include in the Co term, and a freedom 
to shuffle terms between the Di and R away from the exact limit ^ — >■ 0. For our method of dealing with these 
ambiguities, we refer the reader to [35| . which specifies all details relevant for our implementation of the approximate 
NNLO formulas used in the numerical analysis. 

Equations ((!]) and ©-([T]) are written for the case where the factorization and renormalization scales are set 
equal, Hf = /i^. Later on we will consider the case where the running coupling is instead evaluated at an arbitrary 
renormalization scale fXr- To derive such expressions we use 



asifJ-r) 



529^, 
144 



rf 



29 

12^'-^ 



(8) 



with Lrf = ln(/x^//Lij), and re-expand the formulas to NNLO in powers of as(/ir)- 

Although PIM and IPI kinematics are used to describe different differential distributions, one can always integrate 
over the distributions to obtain the total hadronic cross section or the partonic scaling functions fij using the two 
approaches. Since the two kinematics both account for soft gluon emission, they give rise to the same results for fij 
in the limit s — >■ Am^. When s deviates from 4to(, the difference between the two kinematics is formally subleading 



in ^. However, the numerical differences may be visible or even significant if s is much larger than 4toj. In [35[, we 



have shown that the effective-theory predictions for the total cross section from the two kinematics actually agree 
quite well, as long as the exact dependence on the energy of soft gluon radiation is kept in the factorization formula 
([S])-^ The differences between these two kinematics can be regarded as another source of theoretical uncertainty, 
namely that due to power corrections to the soft limit, and used along with scale uncertainties to estimate the total 
uncertainty associated with the calculation of the total cross section at NLO-I-NNLL or approximate NNLO. In the 
following sections we define our procedure for combining the two types of uncertainties and give our final results for 
the total cross section as a function of the top-quark mass in the pole, MS, and IS schemes. 



III. CROSS SECTIONS IN THE POLE SCHEME 



We now specify our method for estimating the theoretical uncertainties from scale variations, PDF variations and 
variations of as in the calculation of the top-pair production cross section, using either PIM or IPI kinematics. We 



^ Previous calculations expanded -Es(^) in the limit z — )■ 1 for PIM kinematics and S4 — > for IPI kinematics, in which case the differences 
between the two types of kinematics become much larger. 



4 



then define a procedure for combining the results obtained in these two schemes. By default, the pole mass mt is used 
in the calculation of the partonic cross sections. We comment later on alternative schemes for defining the top-quark 
mass. 

To estimate the uncertainties associated with scale variations, we view the cross section as a function of the 
renormalization and factorization scales, which by default are chosen as fif — iir — "mt- We then consider two 
methods of scale variations: correlated variations with /^/ = /i^ varied up and down by a factor of two from the 
default value, and independent variations of /x/ and /i,- by factors of two, with the uncertainties added in quadrature. 
We use as our final answer the larger uncertainty from these two methods. 

To combine the results from PIM and IPI kinematics, we first compute the cross sections and scale uncertainties 
in the PIM and IPI schemes separately, and obtain six quantities upiM, ^'''pim' ^'''pim' "'ipIi ^'''ipi' ^''^ipr 
central value and perturbative uncertainties for the combined results are then determined by 

1 r 

Ao-+ = max (ctpim + Actj^jj^ , ctipi + Acr+pj) - a, (9) 
Act" = min (ctpim + AcTpjj^ , cripi + Acj^pj) - a . 

In this way, the central value is the average of the two, and the perturbative uncertainties reflect both the variation 
of the scales and the difference between the two types of kinematics. The PDF uncertainties are estimated as usual 
by evaluating the average of the IPI and PIM results using the PDF error sets at a particular confldence level. 

We quote in Table [J the approximate NNLO predictions obtained with the above procedure at rrit = 173.1 GeV, 
using the 90% CL sets of the MSTW2008 NNLO PDFs To investigate the convergence of the perturbative series, 
we also list the LO and NLO results, obtained using MSTW2008 LO and NLO PDFs, respectively. In the pole scheme, 
the scale uncertainties are generally determined by the correlated scale variations with fir = fJ'f- The one exception 
is the upper error at the Tevatron, which is instead determined by the independent variations of /i^ and /i/ added in 
quadrature. Even though the perturbative uncertainty in the approximate NNLO result includes both scale variations 
and an estimate of power corrections to the soft limit through the difference of IPI and PIM kinematics, it is still 
reduced compared to that in the NLO calculation, which by definition is due only to scale variations. We note that 
the central value and uncertainties of the approximate NNLO results are well contained within the uncertainty range 
predicted by the NLO results, so that the perturbative series to this order is well behaved in the pole scheme. The 
NNLO results are also within the uncertainties of the LO calculation, although the NLO results are slightly higher 
than the LO ones in the case of the LHC. 





Tevatron 


LHC7 


LHC14 


MSTW 


CTEQ 


MSTW 


CTEQ 


MSTW 


CTEQ 


LO 


R cc+2-95+(0.34) 
"•""-1.87-(0.27) 


c .I.+2. 16+0. 33(0. 29) 
^•"^^-l. 42-0. 27(0. 24) 


ir,r,+49+(6) 
-L^^-32-(7) 


-'-'-'^-24-8(7) 


„„,+228+(26) 
"''l-159-(34) 


= =9+157+25(18) 
■^■^^-115-25(19) 


NLO 


„ „+0. 41+0. 47(0. 37) 
^-0.76-0.45(0.24) 


„ __+0. 40+0. 50(0. 43) 
' '-0.74-0.40(0.34) 


,^•0+20+14(8) 
-'-'^^-21-13(9) 


■140 + 18+13(11) 
1^:0-19-12(10) 


000+107+66(31) 
Oo»-106-58(32) 


890+97+41(27) 
°^-96-40(28) 


NNLO approx. 


„ „n+0. 07+0. 63(0. 33) 
"•'-'■^-0.41-0.48(0.25) 


r. q, +0.09+0.53(0.46) 
"•^^-0.44-0.43(0.36) 


1 i:c:+S+14(8) 
^^'-'■^-9-14(9) 


lcq+8+13(ll) 
■3-8-12(10) 


orr + 52+60(30) 
0'-'J-38-59(31) 


0/19+51+40(26) 
''^^-37-40(28) 



TABLE I: Total cross sections in pb for mt = 173.1 GeV with MSTW2008 and CTEQ6.6 PDFs. The first error results from 
the perturbative uncertainty from both scale variations and the difference between PIM and IPI kinematics, the second one 
accounts for the combined PDFs+as uncertainty. The numbers in parenthesis show the PDF uncertainty only. 

For comparison, we also include the results using CTEQ6.6 PDFs 44] in Table HI Since the CTEQ PDFs are based 
on a NLO fit, the same set is used at LO, NLO and approximate NNLO. The statements based on the analysis with 
MSTW PDFs above, including those concerning the moderate size of the NNLO corrections, are also true for the 
analysis with CTEQ PDFs. In this case, however, the LO results at the LHC are significantly lower than the NLO 
and NNLO results. To a certain extent, this shows the potential benefit of switching PDFs as appropriate to the order 
of perturbation theory. On the other hand, LO calculations are usually considered unreliable, so the more important 
observation for the perturbative convergence is the modest size of the NNLO correction. 

The perturbative uncertainties in the approximate NNLO predictions are about the same size at both the Tevatron 
and the LHC. An additional source of uncertainty is related to the experimental value of as{Mz) (where Mz denotes 
the Z-boson mass), which is an input parameter for the running of the strong coupling constant. We estimate this 
uncertainty in combination with the PDF one by employing the method proposed in [45, 46] . Table U shows that the 
uncertainty on as{Mz) adds an error of ±(3 -4)% to the pair-production cross section when the calculation is carried 
out with MSTW2008 PDFs. The error is somewhat smaller, ±(l-2)%, when CTEQ6.6 PDFs are used. The reason 
is that CTEQ6.6 assigns a 90% CL error of ±0.002 to a^^Mz), while for MSTW2008 it is ±0.003. One can conclude 
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that the as{Mz) induced uncertainty is of the same order of magnitude as the perturbative and PDF uncertainties, 
and should not be neglected. 

For an extraction of the top-quark mass through a comparison with the experimental cross section, we also provide 
our results as a function of mt . We parametrize the mass dependence of the approximate NNLO cross section using 
the simple polynomial fit 

a{mt) = Co + cix + C2X^ + c-^x^ + c^x"^ , (10) 

where x = mt/GeW — 173, and Ci are fit coefficients which depend on the collider and the PDF set. The results for the 
fit coefficients including upper and lower errors due to perturbative uncertainties arc shown in Table HIl again using 
MSTW2008 NNLO PDFs. A Mathematica implementation of the fit coefficients can be found with the electronic 
version of this letter, where the combined PDF and uncertainties as well as the fit coefficients using CTEQ6.6 
PDFs are also included. These fits reproduce the approximate NNLO calculations to 1 permille or better in the range 
mt G [150, 180] GeV. For simplicity the uncertainties on fit coefficients are not displayed in Table |lTl When these 
uncertainties are measured in percent of the central value of the cross section, they appear to be roughly independent 
of mt in the range mt € [150, 180] GeV, differing by no more than a percent from those at mt — 173.1 GeV shown in 
Table H 







CO [pb] 


ci [pb] 


C2 [pb] 


C3 [pb] 


C4 [pb] 


Tevatron 


a 

a + Act+ 
a + Act" 


6.64792 X 10° 
6.72257 X 10° 
6.23323 X 10° 


-2.07262 X 10"^ 
-2.09199 X 10"^ 
-1.94555 X 10"^ 


3.61739 X 10"^ 
3.62959 X 10"=^ 
3.40149 X 10"^ 


-4.30451 X 10"^ 
-5.00960 X 10"^ 
-4.03465 X 10-'^ 


8.94347 X lO"'^ 
6.99427 X lO"'^ 
8.17661 X lO"'^ 


LHC7 


a + Aa+ 
a + Aa~ 


1.55546 X 10^ 
1.63325 X 10^ 
1.46870 X 10^ 


-4.66554 X 10° 
-4.92409 X 10° 
-4.40296 X 10° 


8.07632 X 10~^ 
8.50408 X 10"^ 
7.51759 X 10"^ 


-9.93138 X 10"'' 
-1.13686 X 10"^ 
-1.03349 X 10"^ 


1.75303 X 10"^ 
1.66808 X 10"^ 
1.48145 X lO"'^ 


LHC14 


a 

a + Aa+ 
a + Act" 


8.57671 X 10^ 
9.08856 X 10^ 
8.19550 X 10^ 


-2.29929 X 10^ 
-2.44668 X 10^ 
-2.20396 X 10^ 


3.65310 X 10"^ 
3.89564 X 10^^ 
3.54379 X 10^^ 


-4.02700 X 10"^ 
-4.48258 X 10"^ 
-3.81968 X 10^3 


7.41522 X lO"'^ 
8.14287 X 10"^ 
6.25216 X 10"^ 



TABLE II: Fit coefficients in (|10() for the total cross sections with perturbative uncertainties at approximate NNLO, using 
MSTW2008 NNLO PDFs. 



IV. CROSS SECTIONS IN THE MS AND IS SCHEMES 



It is well-known that the pole mass of a quark cannot be defined unambiguously in QCD due to confinement; 
the perturbatively defined pole mass is sensitive to long-distance physics and suffers from renormalon ambiguities of 
order Aqcd [13, El]. In perturbative calculations, the renormalon ambiguity is associated with large higher-order 
corrections to the pole mass, and thus to any observable calculated in this scheme. Therefore, it is worth investigating 
short-distance mass definitions which are free from these shortcomings. In this section, we analyze the cross section 
as a function of the running top-quark mass defined in the MS scheme, and of the threshold top-quark mass defined 
in the IS scheme (37| . 

It is possible to calculate the cross section using the MS mass from the beginning, by performing mass renormal- 
ization in that scheme. However, since we already have the cross section in the pole scheme, it is simpler to convert 
from one scheme to another using t he p erturbative relation between the pole mass and MS mass. This relation is 
currently known to three-loop order [4^. To perform the conversion to the MS scheme, we take that result for QCD 
with five active flavors and write it in the form 



mt = 



1 



as(Mr) 



(11) 



where 



^ 8.23656 



379 



L'rr 



37 

24^ 



23 
12 



(12) 
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with Lra = ln(/x^/m^(/i)) and L,. = ln(/i^//i^). We then decompose the NNLO cross section in the pole scheme as 



0'NNLo(™t) 



as (Mr) 



1 3 



a'^^\Tnt,fJ,r) + 



(13) 



ehminate mt through the relation (jlll) . and re-expand the result in powers of as{fir)- The resulting cross section in 
the MS scheme can be written as 



crNNLo('7l) 

where 



n 3 



cr(2)(m(^),/x,/Xr) , 



(14) 



dmt 



(15) 



mt=m{p) 



m(/i) 



dmt dmt 2 



dmt 



- mt=m(p.) 



The derivatives can be taken either at the level of the hadronic cross section, using fits such as the one in (ITU|) . or 
at the level of the differential cross section before carrying out the phase-space integrations. We have checked our 
calculations by verifying the agreement between the two methods. We note that our method of converting results 
from the pole scheme to the MS scheme is similar to that used in 0, [131 • Indeed, our approximate NNLO results in 
the MS scheme for the choice jl — rn agree with those in the HATHOR program [50| , apart from the piece related to 
the NNLO correction a^^\ which is of course different since we are not working in the s — ^ 4m( limit of the partonic 
cross section. 

Our procedure for combining the results from IPI and PIM kinematics in the MS scheme is analagous to that for 
the pole scheme described above. In the present case, we use by default fj,f = = m(jn). We must also specify the 
scale in the running top-quark mass, for which we use fi = fn.^ We provide results for the cross sections as a function 
of m{m) using the fit 



Co -I- cix + C2X^ -I- c-iX^ -I- cax'^ . 



(16) 



where x = m/GeV - 164. The fit coefficients for the different colhders using the MSTW2008 NNLO PDFs can be 
found in Table Hill those including combined PDF and as uncertainties also with CTEQ6.6 PDFs are included in the 
Mathematica notebook mentioned above. 







Co [pb] 


ci [pb] 


C2 [pb] 


C3 [pb] 


C4 [pb] 


Tevatron 


a 

a + Act+ 
a + Aa~ 


6.66715 X 10° 
6.77748 X 10° 
6.26205 X 10° 


-2.17800 X 10~^ 
-2.21151 X 10"^ 
-2.05259 X 10"^ 


3.95994 X 10"^ 
4.04717 X 10"^ 
3.81108 X 10"^ 


-5.14404 X lO-"^ 
-5.09432 X 10-"' 
-3.88796 X 10-"' 


1.09983 X 10-® 
1.11678 X 10-® 
1.28549 X 10-® 


LHC7 


o 

a + Act+ 
cr -1- Act" 


1.57441 X 10^ 
1.66413 X 10^ 
1.48389 X 10^ 


-4.94191 X 10° 
-5.20036 X 10° 
-4.62721 X 10° 


9.00990 X 10-2 
9.48216 X 10-2 
8.26628 X 10-2 


-9.38583 X 10-'' 
-1.04597 X 10-3 
-1.06324 X 10-3 


2.97762 X lO-'^ 
2.81345 X lO-'^ 
2.22025 X 10-^ 


LHC14 


a 

a + Ao-+ 
a -1- Act" 


8.64542 X 10^ 
9.20794 X 10^ 
8.03504 X 10^ 


-2.42364 X 10^ 
-2.59880 X 10^ 
-2.21450 X 10^ 


3.98093 X 10-^ 
4.34218 X 10-^ 
3.60816 X 10-^ 


-4.89960 X 10-3 
-5.00638 X 10-3 
-4.19956 X 10-3 


8.44709 X lO-'^ 
1.02994 X 10-" 
7.84995 X 10-^ 



TABLE III: Fit coefficients H16|) for the cross section with perturbative uncertainties at approximate NNLO in the MS scheme, 
using MSTW2008 NNLO PDFs. 



2 Variations of p, around values close to m, which would correspond to sampling over different mass definitions, could potentially be used 
as an additional means of estimating systematic uncertainties. However, a numerical analysis shows that our approximate NNLO results 
are very stable for variations of p, around the default value. 
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The results for rn{m) = 164.1 GeV, which corresponds to nit = 173.1 GeV when using the two-loop conversion 
between the pole and MS masses, are shown in Table HVl for MSTW2008 and CTEQ6.6 PDFs. As in the pole scheme, 
we switch the order of the MSTW PDFs according to the order of perturbation theory at which we are working, while 
the CTEQ PDFs are the same in both cases. In the MS scheme, the uncertainties from scale variations are dominated 
by the scheme where fif and /i^ are varied independently, rather than the scheme with correlated = fJ-f variations, 
as was the case in the pole scheme. 





Tevatron 


LHC7 


LHC14 


MSTW 


CTEQ 


MSTW 


CTEQ 


MSTW 


CTEQ 


LO 


o o9+3.91+(0.44) 
°-°^-2.4S-(0.35) 


y r, ,+2.86+0.46(0.40) 
' •^^-1.89-0.38(0.32) 


^gQ+64+(8) 
^""-42-(9) 


, Ol +45+11(9) 
J-'Jl--31-10(8) 


Oyr+291+(32) 

°'>^-204-(43) 


y(^r + 199+30{21) 
'"0-145-30(23) 


NLO 


y oq+0-ll+0-50(0.40) 
' ••^'^-0.49-0.47(0.25) 


y OQ+0. 10+0.57(0.50) 
' ••3^-0.48-0.45(0.39) 


1 70+11 + 15(10) 
1 ' ''-19-14(10) 


1R7+10+15(12) 
1"' -17-13(11) 


QQ, +79+71(35) 
»»J-_96-62(36) 


Q9C+71+44(29) 
»^J_87-43(31) 


NNLO approx. 


f, „,+0. 11+0.58(0.33) 
"•"^-0.40-0.43(0.23) 


fi q9+0-12+0.52(0.46) 
"••'^-0.43-0.42(0.37) 


1 c^7+9+13(8) 
1^' -9-13(9) 


1 c:zl+^'+13(ll) 
l'^^-9-12(10) 


Of. 9+56+54(30) 
00^-61-53(32) 


0^0+56+37(26) 
0^0-61-38(28) 



TABLE IV: Total cross sections in pb in the MS scheme, for m(m) = 164.1 GeV. The first error results from the perturbative 
uncertainty from both scale variations and the diflference between PIM and IPI kinematics, the second one accounts for the 
combined PDFs+as uncertainty. The numbers in parenthesis show the PDF uncertainty only. 

We observe that the results obtained from the approximate NNLO formulas are quite close to those in the pole 
scheme shown in Table HI both in the central values and in the errors. Given this good agreement, which is roughly 
independent of the exact value of the top-quark mass as shown by the fits, it makes little practical difference whether 
one extracts the pole mass using the approximate NNLO results, and then determines the MS mass using the per- 
turbative conversion (jlip. or whether one determines the MS mass directly, using the experimental results along with 
the fits at approximate NNLO. This statement would not be true at very high orders in perturbation theory, since the 
renormalon ambiguity inherent to the pole mass would lead to large corrections not present in a short-distance scheme 
such as the MS scheme. But given the present accuracy of perturbative calculations and experimental measurements, 
this does not yet appear to be an issue. 

It is of course still interesting to study whether even at low orders the perturbative expansion is better behaved 
in the MS scheme than in the pole scheme. We observe that the perturbative uncertainties at NLO are generally 
smaller in the MS scheme than in the pole scheme, and that the central values are relatively higher compared to the 
approximate NNLO calculation. For this reason, the overlap between the NLO and approximate NNLO results is 
actually better in the pole scheme than in the MS scheme."^ These results differ from those obtained in the s — >■ Am^ 
limit, where the approximated NNLO corrections and the perturbative uncertainties at that order are significantly 
smaller in the MS scheme than in the pole scheme Q. 

To elaborate further on these results, we note that the re-organization of the perturbative expansion in the MS 
scheme compared to the pole scheme is accomplished by the terms in square brackets in (|15l) . To understand whether 
these terms are expected to cancel againt unphysically large corrections in the pole scheme, we note that the main 
source of mass dependence in the Born level cross section is due to phase-space factors: the lower limit of integration in 
^ , and an overall factor of -^/l — Amf/s in the partonic cross section related to two-body phase space and multiplying 
the Born-level matrix element. The derivatives contained in the terms in square brackets are mainly sensitive to those 
sources of mt dependence. However, the phase space of the pair production is more indicative of the pole mass 
than of an MS mass. Indeed, we are calculating the cross section for on-shell quarks according to the narrow width 
approximation. If the cross section is instead calculated in the MS scheme, the terms in the square brackets of the 
NLO and NNLO pieces of ([15]) give sizeable negative corrections, which are accounted in the pole scheme by using a 
numerically higher value of the mass in the LO and NLO cross sections. Since the most appropriate mass scheme for 
a given process is the one where the higher-order corrections are expected to be smallest on physical grounds, it does 
not seem to us that the MS scheme is the optimal choice for this case. 

As an alternative to the MS mass, we consider the group of short-distance masses known as threshold masses 
[36j . At lower orders in perturbation theory, these are closer numerically to the pole mass, but they do not suffer 
from renormalon ambiguities at higher orders. The cross section in these schemes can be easily calculated from the 
pole-scheme results, using an analogous procedure to the MS scheme calculation. It is evident that at approximate 
NNLO the numerical difference between these results and the MS and pole-scheme results will be quite small once 
the numerical value of the mass is adjusted appropriately, but we nonetheless illustrate this with a specific example. 



3 The overlap between LO and NLO is worse at the Tevatron and improved at the LHC compared to the pole scheme, but as mentioned 
earlier we consider the more important issue the overlap between the NLO and approximate NNLO results. 
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In particular, we consider the cross section as a function of the IS mass introduced in [37[. The IS mass is defined 
through the perturbative contribution to the mass of a hypothetical n = 1, '^Si to pon ium bound state. To perform 
the conversion to this scheme, we write its relation with the pole mass in the form 



rut 



as 



TT 9 V TT 



2 , ,/23, ?,^lr 

9 \ d 4as(/ir)w 



181 

Is" 



9 




and follow the same procedure as for the MS scheme calculation with the appropriate replacements, cf. ([Tl]) . Note 
that in the above relation nag is counted as 0(1) and is not expanded. The results are listed in Table |V] for the value 
= 172.3 GeV, which corresponds to a pole mass of mt = 173.1 GeV using the two-loop conversion above. The 
approximate NNLO results in this scheme are very similar to those in the pole and MS schemes, but the moderate size 
of the NNLO correction is more indicative of the pole scheme than of the MS scheme. This leads us to conclude once 
again that although at yet higher orders in perturbation theory the pole mass would be disfavored, at approximate 
NNLO accuracy this is not yet a problem. 





Tevatron 


LHC7 


LHC14 


MSTW 


CTEQ 


MSTW 


CTEQ 


MSTW 


CTEQ 


LO 


f. oq+3-02+(0.35) 
O-O'J_i.92_(0.28) 


r PQ+2. 21+0. 34(0. 30) 
O-0y_l. 46-0.28(0.24) 


,r,.+50 + (6) 

-^^*-33-(7) 


,„o+35+9(7) 
-^^■^-24-8(7) 


„Q„+223+(26) 
"^"-163-(34) 


rp. + 160+25{18) 
>J"'±_117_25(19) 


NLO 


f- 09+0'39+0.48(0.38) 
°-°^-0. 75-0. 46(0. 24) 


r. o„+0. 38+0. 51(0. 44) 
"•°'-0. 73-0. 41(0. 34) 


if;9+19+14(9) 
-^"^-21-13(9) 


1 t;n+"+14(ll) 
-^■-"^-19-12(10) 


q|-, 9+106+66(32) 
^^^-106-59(33) 


o,, +96+41(27) 
0^-^-96-40(29) 


NNLO approx. 


r. f,r+0. 06+0.63(0.32) 
"•""^-0. 38-0. 47(0. 24) 


r. qo+0. 08+0. 54(0. 47) 
^•^■^-0.40-0.42(0.36) 


lc;c+7+14(9) 
-^"-'"_8-14(8) 


lc4+7+13(ll) 
-^■^^-8-12(10) 


orq+47+59(30) 
°>J'^_35_58(32) 


04C+46+39(25) 
°^"-35-40(29) 



TABLE V: Total cross sections in pb in the IS scheme, for mt = 172.3 GeV. The first error results from the perturbative 
uncertainty from both scale variations and the difference between PIM and IPI kinematics, the second one accounts for the 
combined PDFs+Qs uncertainty. The numbers in parenthesis show the PDF uncertainty only. 



V. CONCLUSIONS 



We have presented predictions for the total inclusive cross section for top-quark pairs at hadron colliders at ap- 
proximate NNLO in QCD. Our calculations are based on soft gluon resummation to NNLL order in PIM and IPI 
kinematics, carried out within the context of effective field theory. They represent the state-of-the-art, combining all 
knowledge presently available about higher-order QCD corrections to the production cross section. The perturbative 
uncertainties associated with our results are estimated in two ways: through the standard method of variations of 
factorization and renormalization scales, and also through the difference between the two types of kinematics. The 
latter gives a means of estimating the size of perturbative power corrections to the soft limits in which the approximate 
NNLO formulas are derived. The results presented here consolidate those previously presented in [H, |33 • We have 
also provided a computer program which calculates the total cross section within our approach. 

The total production cross section can be used along with experimental measurements to extract the top-quark 
mass. An advantage of such extractions is that the theory calculations are carried out in a well-defined renormalization 
scheme for rat- For very precise extractions of the top-quark mass the pole mass is disfavored, because it is only defined 
up to a renormalon ambiguity of order Aqcd- In practice, however, we have not observed a poor convergence of the 
perturbative series up to NNLO in the pole scheme compared to the MS scheme, and pointed out that the group of 
short-distance masses known as threshold masses may be equally appropriate. We have provided numerical fits of our 
results as a function of the mass in both the pole and MS schemes, including perturbative and PDF uncertainties, in 
addition to those from the strong coupling constant, which are non-negligible at this level of accuracy. 

The results presented in this letter can be used directly by the experimental collaborations at the Tevatron and 
LHC in top-quark mass measurements from the corresponding production cross sections. 
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